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We report the observation of the B meson decay B ± — > J /ipnK ± and evidence for the decay 
B° -f J/ipr}K%, using 90 million BB events collected at the T(AS) resonance with the BABAR 
detector at the PEP-II e + e~~ asymmetric-energy storage ring. We obtain branching fractions of 
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B(B ± -> J/i/;r 7 ii' ± )=(10.8±2.3(stat.)±2.4(syst.)) X 1(T 5 and B(B° -> J/>r7ifg)=(8.4±2.6(stat.)± 
2.7(syst.)) x 10~ . We search for the new narrow mass state, the X(3872), recently reported by 
the Belle Collaboration, in the decay B ± — > X(3872)K ± , X(3872) — » J/tpr] and determine an upper 
limit of B(B ± -> X(3S72)K ± -> J/ipvK*) < 7.7 x 10" 6 at 90% C.L. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The study of charmonium states produced in exclusive 
B meson decays led to observations of known charmo- 
nium states and recently to the discoveries of new states. 
Since B mesons can decay via color-suppressed b — ► ccs 
quark transitions, the charmonium states are typically 
produced in final states with kaons. Many known char- 
monium states have been observed in decays such as B — > 
J/if)K(*\ ip(2S)K(*\ XcK^, an d Vc (lS)K^ and evi- 
dence for new states such as a candidate for the 7y c (3654) 
has been published Recently the Belle Collabora- 
tion observed a new narrow mass state with a 3.872 
GeV/c 2 mass produced in the decay B ± -> X (3872) K ± , 
X (3872) — > 7r + 7r _ J/ip. This new state may be the hith- 
erto undetected J PC = 2~~ 1 3 I?2 charmonium state Q. 
However, such a state should have a large radiative El 
dipole transition into 7Xd, which Belle does not observe, 
and theoretical models |3j predict a smaller mass split- 
ting, relative to the 7^(3770), than observed. Unconven- 
tional explanations include a molecule |4j formed with 
charmed D and D* mesons, since the X(3872) has a 
mass exactly at D*° (2007) + D° (1864) threshold. Al- 
ternatively, this new state may be a hybrid charmonium 
state [5j formed of cc ± gluons since color octet charmo- 
nium states may be produced in exclusive B decays Q. 

To further elucidate the nature of the X(3872), we 
performed an analysis on the new exclusive decay B — » 
J/^K, to search for X (3872) -> J/ipr). If the X (3872) is 
a conventional charmonium state, its decays may be sim- 
ilar to the ip(2S), which decays into J/ipir + ir~ and, with 
a factor ten smaller relative rate, into J/ipT). If instead, 
it is a hybrid charmonium state, it is also predicted || to 
decay into J/ipinr and J/tpr] with possibly an enhanced 
rate in the rj channel. 

The decay B — * J/ipnK is similar at the quark level 
to other color-suppressed decays such as B — > J/ipcpK 
which has been observed with a branching fraction of 
(4.4 ± 1.4 ±0.5) x 10~ 5 0. Hence it might be expected 
that B — > J/tpr]K has a comparable branching fraction. 

The data used in this analysis correspond to a total 
integrated luminosity of 81.9 fb _1 taken on the T(4S) 
resonance, producing a sample of 90.0 ±1.0 million 
BB events (N B g). Data were collected at the PEP-II 
asymmetric-energy e + e~ storage ring with the BABAR 
detector, fully described elsewhere [8j. The BABAR de- 
tector includes a silicon vertex tracker and a drift cham- 
ber in a 1.5-T solenoidal magnetic field to detect charged 
particles and measure their momenta and energy loss. 
Photons, electrons, and neutral hadrons are detected in 
a CsI(Tl)-crystal electromagnetic calorimeter. An inter- 



nally reflecting ring-imaging Cherenkov detector is used 
for particle identification. Penetrating muons and neu- 
tral hadrons are identified by resistive-plate chambers in 
the steel flux return. Preliminary track-selection crite- 
ria in this analysis follow previous BABAR analyses 
and the detailed explanation of the particle identification 
(PID) is given elsewhere 0, [Icj . 

The intermediate states in the charged (J/tpr/K^) and 
neutral (J/iprjKg) modes used in this analysis, J/ip — > 
e + e~, J ftp — > r\ — » 77 and Kg — > 7r + 7r~, are se- 

lected within the mass intervals 2.95 < M(e + e~) < 3.14, 
3.06 < M(n + n~) < 3.14, 0.525 < M( 77 ) < 0.571, and 
0.489 < M (tt+tt-) < 0.507 GeV/c 2 . The mass interval 
for e + e _ is larger than than for to enable detection 

of events with Bremsstrahlung in the detector. The K s 
decay length in the lab frame is required to be greater 
than 0.1 cm. 

Determination of the signal and the background uti- 
lizes two kinematic variables 7]: the energy difference 
AE between the energy of the B candidate and the 
beam energy E£ in the T(iS) rest frame; and the beam- 
energy-substituted mass toes = \J (E^) 2 — (p^) 2 , where 
is the reconstructed momentum of the B candidate 
in the T(45) frame. Signal events should be concen- 
trated in a rectangular signal-box region bounded by 
b^ES ~ m B| < 7.5 MeV/c 2 , where me is the mass of 
B meson and \AE\ < 40 MeV. 

Before the data were analyzed, the final selection cri- 
teria were optimized separately for the charged and neu- 
tral modes using a Monte Carlo (MC) simulation of the 
signal and the known backgrounds. Motivated by the 
B — > J/ipcpK measurement, the ab initio value of the 
branching fraction for B — > J/iprjK used in the signal 
MC was 5 x 10 -5 . The number of reconstructed MC 
signal events n™ c and the number of reconstructed MC 
background events n™° in the signal-box were used to 
estimate the sensitivity ratio, n™ c / \Jn™ c ± n™ c . This 
ratio was maximized by varying the selection criteria on 
the rj mass, a 7r° veto, the photon helicity angle from 
the rj decay and the thrust angle. The 77 mass inter- 
val of the 77 candidate as specified earlier was chosen by 
this procedure. In the charged (neutral) mode, if either of 
the photons associated with an 77 candidate, in combina- 
tion with any other photon in the event, forms a 77 mass 
within 17(10) MeV/c 2 of the nominal 7T° mass, the rj can- 
didate is vetoed as a ir° background. The 77 candidate is 
rejected if |cos^| is greater than 0.93(0.81), where 9^ is 
the photon helicity angle 9] in the rj rest frame. Signal 
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events have a uniform cos 9^ distribution whereas combi- 
natorial background of random pairs of photons typically 
has a distribution that peaks near ±1. 

To separate two-jet continuum events from the more 
spherical decays of B mesons produced nearly at rest 
from T(AS) — > BB, the angle 9t between the thrust 
direction of the B meson candidate and the thrust di- 
rection of the remaining charged tracks and photons in 
the event is calculated. We reject events when |cos#t| is 
greater than 0.8(0.9), since the distribution in cos#t is 
flat for BB events, while background e + e~ — + qq contin- 
uum events peak at cos#t = ±1. 

The data, after these cuts, are shown in Figs. and [21 
where (a) is a scatter plot of AE versus to.es, (b) is the 
AE histogram and (c) is the toes histogram (solid line). 
We find evidence for B signals in both the J/iprjK^ and 
J/ipTjKg modes. 
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FIG. 1: For B — > J/ipriK ± , the AE versus 7ties event distri- 
bution (a) is shown with vertical and horizontal bands defined 
by limits, |m E s -m B | < 7.5 MeV/c 2 and \AE\ < 40 MeV, re- 
spectively. The intersection of these bands corresponds to the 
signal-box region defined in the text. The AE projection (b) 
is shown for events in the vertical band that contains the tues 
signal region. The tties projection (c) is shown for events in 
the horizontal band that contains the AE signal region. The 
dashed histogram represents the estimated background and is 
described in the text. 



To determine the branching fraction for these modes, 
we first find the number of signal events, which is defined 
as n s = riQ — rtb, where no is the number of events in 
the signal-box region, and «b is the estimated number of 
background events. For each mode, rib is obtained from 
fitting the toes distribution for events with \AE\ < 40 
MeV with the line shape of a Gaussian function and an 
ARGUS function || , which is an empirical parameteriza- 
tion of the background shape. The fit parameters are the 
normalization and mean of the Gaussian and the normal- 
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FIG. 2: The AE and m E s distributions for B° J/ij)r]K%. 
The descriptions of Figs. 2(a), (b), and (c) follow those of 
Figs. 1(a), (b), and (c), respectively. 



ization of the background curve. The width of the Gaus- 
sian is fixed to the value determined by MC simulation 
and the shape of the background curve is fixed to a best 
fit to the data toes distribution with the AE sideband 
region of 0.10 < \AE\ < 0.14 GeV for the B ± mode and 
0.08 < \AE\ < 0.28 GeV for the B° mode. Figs.OJc) and 
|2Ic) show the resulting Gaussian and background curves 
(solid) and the background events (dashed histogram) 
from the AE sideband regions normalized to the data in 
the signal region. Integrating the background curve over 
the signal-box region we obtain rib and its uncertainty, 
(7b. Results are listed in Table [i] Additional checks on 
the background shapes were performed. Using data, a 77 
mass sideband that is outside the nominal 77 mass was se- 
lected and a similar toes background shape was found. 
Using MC simulations of inclusive B J /if) backgrounds 
another toes background shape was obtained and fit. If 
we include a fit component due to this background, which 
is well described by a broad Gaussian distribution, and 
we refit the data sideband distribution, the background 
results do not change. 

The branching fraction is calculated as B — n s /(N B gX 
ex/) where e is the efficiency and / is the product 
of secondary branching fractions for the J/ip, rj, and 
Kg. Efficiencies are determined by MC simulation with 
three-body phase space and the branching fractions of 
T(45) -> B+B- and T(AS) -> B°B° are assumed to be 
equal. Results on B are given in the last column of Ta- 
ble [I] where the first and second errors are statistical and 
systematic, respectively. The statistical error is derived 



from the uncertainty in n s which is -\/ n o + o^. 

The systematic error, cr sys , for each mode (charged 
/neutral) is determined by adding in quadrature the 
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TABLE I: Efficiencies, number of signal-box and background events, 90% C.L. of the number of events and the branching 
fraction upper limits, P-values and branching fractions 



Mode 


e 


n n b ± cr b 


^90% 


90% C.L.U.L. 


P-value 


Branching Fraction 




10.75% 


99 50.3 ±3.0 


70.0 


< 15.5 x 10~ 5 


2 x 10~ 8 


(10.8 ±2.3 ±2.4) X 10" 5 


J/iln)K% 


8.53% 


39 18.5 ± 1.7 


34.5 


< 14.1 x 10~ 5 


9 x 10~ 5 


(8.4 ± 2.6 ± 2.7) x 10" 5 



percentage uncertainty on each of the following quan- 
tities: N B g (1.1/1.1); secondary branching fractions 
(2.48/2.52); MC statistics (1.77/2.17); PID, tracking, 
and photon detection efficiencies (8.2/8.3); ir° veto 
(8.1/8.3); r) mass range (3.40/3.14); background param- 
eterization (16.7/27.0); and model dependence (5.1/9.5). 
The total systematic errors for the charged and neutral 
modes are 22.0% and 32.0%, respectively. The uncer- 
tainties in the PID, tracking, and photon detection ef- 
ficiencies are based on the study of data control sam- 
ples 0- The uncertainty in the 7r° veto efficiency was 
studied by measuring the veto efficiency on the inclusive 
r\ rate in data and MC. The uncertainty due to the r\ mass 
selection was determined by comparing the measured rj 
mass resolution in inclusive 77 decays to the 77 mass resolu- 
tion from the signal MC. The background parameteriza- 
tion uncertainty was estimated by changing the ARGUS 
shape parameter by ±1 standard deviation, refitting the 
771ES data distribution, and recalculating the number of 
signal events. Although this analysis used MC events 
generated with three-body phase space to determine the 
final efficiencies, additional systematic uncertainties due 
to the decay model dependence are estimated. The effi- 
ciency uncertainty due to unknown angular distributions 
and intermediate resonances has been estimated by com- 
paring the efficiencies obtained in five different MC gener- 
ated models. These include 100% transversely polarized 
J/ip, 100% longitudinally polarized J/ip, large two-body 
J/iprj mass, large two-body rjK mass and small two-body 
J/ipK mass. The resulting relative change in efficiencies 
was used to estimate the production model uncertainty. 
The resulting total a sys for each mode is used to deter- 
mine the B systematic errors in Table [I] 

The P-value for null hypothesis (no signal) is the Pois- 
son probability that the background events fluctuate to 
> no. Assuming the probability distribution function of 
the background is a Gaussian with mean and stan- 
dard deviation Ob, we calculate the Poisson probabilities 
with different background values weighted by this Gaus- 
sian distribution to determine the final P-value for each 
mode. The resulting P-values are equivalent to a statis- 
tical significance of 5.6cr and 3.9cr for the charged and 
neutral modes, respectively. 

We also determine the 90% confidence level upper limit 
(C.L.U.L.) on the branching fraction using no, rib, and 
<7b, in the signal region, and <7 sys . The Bayesian upper 
limit on the number of signal events, Ng %, is obtained 



by folding the Poisson distribution with two Gaussian 
distributions representing the background and system- 
atic uncertainties and integrating the resulting function 
to the 90% confidence level (C.L.). This assumes that 
the a priori branching fraction distributions are uniform. 
The charged and neutral results, J/ipirK^ and J/tprjKg, 
are listed in Table [I] 

Our resulting branching fractions are comparable to 
the color-suppressed decay B — > J/ip(j)K branching frac- 
tion. The ratio of the charged (J/tpriK^) to neutral 
(J/ipijK^) branching fractions is consistent within errors 
to the expected value of two. 

We search for the X(3872) in B -> XK, X -> J/ipr) 
now only selecting the signal region, mEs — rn&\ < 
7.5 MeV/c 2 and |AS| < 40 MeV. The resulting J/i/jt] 
mass distribution is shown in Figure |3| The two-body 
mass resolution from Monte Carlo studies is 6 MeV/c 2 . 
There is evidence for the ip(2S) and no evidence for 
the X(3872). Using the measured branching fraction 
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FIG. 3: The J/ipV mass distributions from B — ► J/ipr]K 
and B° -> J/ipr/Ki. The arrows indicate where the ip(2S) 
and X(3872) signals would appear. 



B(B ± -» ^S)^) = (6.8 ± 0.4) x 10" 4 H3, we ex- 
pect to reconstruct 12 ± 1 events in the charged mode in 
the J/tpri mass region below 3.710 GeV/c 2 and we observe 
15. After restricting the mass to 3.85 < M(J/iprj) < 3.89 
GeV/c 2 , we fit the tbes plot with the same procedure 
as before and obtain an upper limit for the product 
branching fraction B{B ± -> X(3872)K ± , X -> J/tpr)) 
< 7.7 x 10~ 6 at 90% C.L. 

Our resulting upper limit may be compared to the 

Rpllp rpcnlt II B(B ± ^X(3872)K ± ^J/^ + 7T-K ± ) _ , f , „ , 

Uelle result y, e(B ±^ l/ , (2S )K±->j/,/>?r+ 7 r-K±) ~ i b - d ± 
1.2 ± 0.7)%. Using B{B ± -> •tM2S)K ± J/^tt+tt- K±) 
= (2.0 ± 0.15 ± 0.22) xl0~ 4 11] it can be deduced that 
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B(B ± -> X(3872)j£: ± -> J/ipw+n-K^ (12.6 ± 2.8 ± 
1.2) x 1(T 6 . If the matrix elements for AT(3872) -> 
J/ipir + TT~ and J/V"7 are similar to those of the ip(2S) 
and we include the larger phase space for the decay of 
X(3872) — » J/V"7 relative to the ip(2S), then we would 
expect S(B ± -> AT(3872)iC b -► J/^rjK^ - 3 x 10~ 6 . 
Our upper limit is within a factor two of this estimate. 
This result is not in contradiction with the charmonium 
interpretation of the X(3872). 

In conclusion, we observe the new decay mode B — > 
J/iprjK with branching fractions of B(B ± — > J/ipr/K^) = 
(10.8 ± 2.3 ± 2.4) x 1CT 5 and B(B° -> J/ipr]K%) = (8.4± 
2.6±2.7) x 10~ 5 . We set an upper limit for the A(3872) in 
the product branching fraction, B{B ± -> X(3872)K ± -> 
J/ipr/K*) < 7.7 x 1(T 6 at 90% C.L. 
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